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(1)  Phase  retrieval  by  projections:  In  [10]  we  see  the  solution  to  a 
50  year  old  engineering  problem  with  applications  to  x-ray  crystal¬ 
lography,  x-ray  tomography,  crystal  twinning,  and  much  more.  The 
results  here  on  phase  retrieval  by  projections  are  truly  significant. 
Previous  work  indicated  that  it  would  take  nlog(n)  arbitrary  rank 
projections  to  do  phase  retrieval.  We  showed  that  2n-l  projections 
is  sufficient  in  the  real  case  and  surprised  everyone  even  more  by  giv¬ 
ing  evidence  that  actually  fewer  than  2n-l  projections  might  work 
-  an  unimaginable  result  with  the  prospect  of  major  applications. 
We  are  continuing  work  in  this  direction  to  find  the  least  number  of 
projections  needed  and  to  do  the  complex  case  and  are  just  starting 
to  work  directly  with  the  engineers  for  implementation. 

(2)  In  [11]  we  introduced  yet  another  new  area  of  research  in  frame  the¬ 
ory:  Integer  frames.  The  idea  here  is  to  speed  up  calculations  and 
accuracy  in  applications  of  frame  theory.  Integer  frames  are  Hilbert 
space  frames  for  which  all  their  coefficients  are  integers.  The  advan¬ 
tages  here  include:  (a)  We  do  not  need  to  approximate  the  frame 
coefficients  speeding  up  calculations,  (b)  Quantization  sometimes 
becomes  a  problem  when  there  are  multiple  levels  of  processing  and 
this  eliminates  one  level  of  quantizations  for  the  process.  This  is 
an  ongoing  project  because  many  of  the  examples  constructed  in 
[11]  are  sparse  frames  and  we  really  need  non-sparse  frames  -  and 
preferable  something  close  to  equiangular  frames. 

(3)  In  [12]  we  made  the  first  systematic  study  of  outer  product  frames. 
These  frames  have  the  potential  for  serious  applications  but  are  very 
difficult  to  work  with.  One  major  result  shows  that  equiangular 
tight  frames  give  the  absolute  best  frame  bounds  for  outer  product 
frames.  We  expect  this  paper  to  open  up  this  area  to  serious  ap¬ 
plications.  This  is  also  an  ongoing  project.  There  are  a  number  of 
critical  construction  problems  left  here  which  need  to  be  resolved  in 
order  to  make  this  class  of  frames  really  usable. 

(4)  In  [9],  we  made  the  first  systematic  study  of  the  distribution  of 
frame  coefficients  of  Hilbert  space  frames.  This  is  fundamental 
to  this  field  and  should  have  been  done  20  years  ago,  but  no  one 
had  successfully  tackled  this  area.  Our  work  is  comprehensive  and 
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invaluable.  This  project  is  also  ongoing  because  we  need  to  do  ran¬ 
dom  methods  here  where  the  outcomes  will  be  significantly  better 
than  currently  available. 

(5)  In  [2]  we  surprised  everyone  by  showing  that  frame  expansions 
of  Parseval  frames  are  1-unconditional  series  and  then  giving  a  com¬ 
plete  classification  of  frames  with  this  property.  This  project  is  com¬ 
plete. 

(6)  In  [3]  we  introduce  a  new  topic:  weaving  frames.  This  topic  has  the 
potential  for  important  applications  in  pre-processing  signals  and 
dealing  with  multiple  wireless  sensor  networks.  This  paper  is  a  com¬ 
prehensive  study  and  should  generate  parallel  research  and  applica¬ 
tions.  This  is  an  ongoing  research  project  since  for  applications  we 
need  to  understand  weaving  for  Gabor  frames. 

(7)  [4]  is  a  continuation  of  the  paper  by  Bodmann  Cahill  and  Casazza, 
Fusion  frames  and  the  restricted  isometry  property.  (English)  Nu- 
mer.  Funct.  Anal.  Optim.  33,  No.  7-9,  770-790,  2012. 

In  analogy  with  this  preceding  paper,  the  correspondence  between 
RIP  matrices  and  fusion  frames  obtained  by  a  partial  orthonormal¬ 
ization  strategy  is  at  the  core  of  the  results.  In  contrast  to  the  preced¬ 
ing  paper,  a  specific  construction  of  RIP  matrices  by  random  Gauss¬ 
ian  matrices  is  used  in  order  to  derive  stronger  consequences  for 
the  resulting  fusion  frames  consisting  of  independent,  uniformly  dis¬ 
tributed  subspaces.  A  measure  concentration  argument  shows  that 
the  Hilbert-Schmidt  inner  products  between  the  orthogonal  projec¬ 
tions  onto  the  random  subspaces  concentrate  near  an  average  value. 
Overwhelming  success  probability  for  near  tightness  and  equiangu- 
larity  is  guaranteed  if  the  dimension  of  the  subspaces  is  sufficiently 
small  compared  to  that  of  the  Hilbert  space  and  if  the  dimension 
of  the  Hilbert  space  is  small  compared  to  the  sum  of  all  subspace 
dimensions. 

(8)  The  paper  [5]  paper  leverages  the  results  on  the  nearly  equal  norm 
and  nearly  tight  fusion  frames  to  derive  error  bounds  for  packet 
encodings  in  the  presence  of  data  loss  (erasures).  The  fusion  frames 
encode  a  vector  in  a  Hilbert  space  in  terms  of  its  components  in 
subspaces,  which  can  be  identified  with  packets  of  linear  coefficients. 
The  fusion  frame  performance  is  evaluated  under  the  assumption 
that  the  vector  to  be  transmitted  is  uniformly  distributed  on  the 
sphere  and  when  part  of  the  packets  is  transmitted  perfectly  and 
another  part  is  lost  in  an  adversarial,  deterministic  manner.  The 
performance  is  measured  by  the  mean-squared  Euclidean  norm  of 
the  reconstruction  error  when  averaged  over  the  transmission  of  all 
unit  vectors.  The  main  result  is  that  a  random  selection  of  fusion 
frames  performs  nearly  as  well  as  previously  known  optimal  bounds 
for  the  error,  characterized  by  optimal  packings  of  subspaces,  which 
are  known  not  to  exist  in  all  dimensions. 
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[7]  Shearlets  are  a  type  of  frames  that  provide  directional  multiscale 
representations  with  optimal  sparsity  properties  for  a  class  of  piece- 
wise  smooth  signals  with  piecewise  smooth  singularities.  The  goal  of 
this  paper  is  to  show  that  a  variant  of  shearlets  can  be  constructed 
with  the  help  of  standard  wavelet  filters  and  standard  Gabor  win¬ 
dows.  Gabor  shearlets,  however,  are  based  on  a  different  group  rep¬ 
resentation  than  previous  shearlet  constructions.  Unlike  the  usual 
shearlets,  the  new  construction  can  achieve  a  redundancy  as  close  to 
one  as  desired.  In  combination  with  Meyer  filters,  the  cone-adapted 
Gabor  shearlets  constitute  a  tight  frame  and  provide  low-redundancy 
sparse  approximations. 

[16]  gives  constructions  of  tight  and  full  spark  Chebyshev  frames 
from  truncations  of  Vandermonde-like  matrices  of  orthogonal  poly¬ 
nomials  are  presented.  These  are  frames  with  real  entries.  Worst 
case  coherence  analyses  are  also  carried  out.  We  show  that  for  suffi¬ 
ciently  high  degree,  the  minimum  angle  between  distinct  frame  vec¬ 
tors  is  bounded  below  by  about  44  degrees.  In  a  related  study,  we 
also  provide  a  worst  case  coherence  analysis  for  equal  norm  tight 
M  x  N  frames  from  truncated  DFT  matrix.  The  cosine  of  the  small¬ 
est  angle  between  these  distinct  frame  elements  is  asymptotic  to 
smiaQ')  where  a  =  M/N. 

[17]  In  this  article,  a  fast  and  super  resolution  method  for  direction 
of  arrival  (DOA)  estimations  is  proposed  under  low  signal-to-noise 
ratio  using  a  limited  number  of  snapshots  (of  measurements).  The 
method  is  based  on  a  sparse  signal  reconstruction  technique  of  null 
space  tuning  in  the  context  of  compressed  sensing  and  sparse  repre¬ 
sentations.  A  crucial  correlation  operation  is  proposed  to  mitigate 
the  noise  effect  from  the  sparse  representation  point  of  view.  The 
proposed  method  has  the  characteristics  of  simultaneous  high  res¬ 
olution,  robustness  to  noise,  and  effective  at  estimating  number  of 
sources.  The  algorithm  is  also  computationally  efficient. 

[18]  In  array  signal  processing,  signal  detections  and  direction-of- 
arrival  (DOA)  estimations,  for  many  years,  suffer  from  poor  resolu¬ 
tions  when  several  signals/targets  are  close  to  each  others.  In  these 
problems,  observed  signals  are  naturally  sparse  linear  combinations 
of  array  direction  vectors  sampled  at  arriving  angles.  The  matrix 
A  of  array  direction  vectors  sampled  at  (typically)  large  number  of 
angle  points  is  a  large  redundant  frame  with  columns  being  frame 
vectors.  If  the  dimension  of  the  ambient  vector  space  is  n,  and 
the  column  of  A  is  N.  We  have  shown  that  there  is  a  common 
sparsity-inducing  dual  frame  for  a  set  of  n  linearly  independent  vec¬ 
tors.  Array  direction  matrix  A  is  typically  full  spark.  Consequently, 
sparsity-inducing  duals  can  be  evaluated  a  priori  for  any  blocks  of  n 
columns  of  A.  It  becomes  therefore  particularly  effective  for  DOA 
estimations  when  signals  are  clunked  together.  This  approach  does 
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not  need  the  typical  t \  -ruin  procedure,  as  applying  sparse  duals  to 
observed  signals  yields  directly  the  few  non-zero  coefficients  (signals) , 
under  a  multiresolution  pursuing  strategy. 

[19]  This  article  is  about  a  group  of  fast  iterative  and  threshold¬ 
ing  algorithms  for  sparse  signal  recovery.  The  algorithms  are  faster 
and  more  effective  than  most,  if  not  all,  known  algorithms.  Ana¬ 
lytical  convergence  results  are  also  established.  One  core  algorithm 
is  shown  to  converge  in  finite  many  steps.  Quite  convincing  exam¬ 
ples  are  shown  through  numerical  experiments.  The  algorithms  are 
particularly  effective  for  large  scale  problems. 

[20]  This  article  introduces  and  studies  the  notion  of  sparse  dual 
frames  for  the  first  time.  Sparse  dual  frames  in  this  context  refer 
to  dual  frames  of  smaller  or  smallest  nonzero  vectors  and  nonzero 
entries  within  a  vector.  Theoretical  lower  bound  of  the  sparsity  for 
random  frames  and  Gabor  dual  frames  are  obtained.  Gabor  dual 
frames  of  the  smallest  time  and  frequency  support  are  obtained. 
Analysis  of  sparse  dual  Gabor  frames  results  in  new  duality  results 
which  are  monumental  for  the  sparse  dual  frame  analysis  in  the  case 
of  dual  Gabor  frames. 

[21] 

When  signals  have  sparse  (coherent)  frame  representations,  i.e. , 
signals  are  sparse  with  respect  to  coherent  frames,  the  compressed 
sensing  problem  becomes  much  more  complicated.  We  have  seen 
that  a  sparse-dual-frame  based  .^-analysis  approach  for  the  sparse 
signal  recovery  is  the  most  effective  one  among  all  known  methods, 
such  as  the  l\  synthesis  approach  and  the  conventional  G-analysis 
approach.  Meantime,  we  also  have  an  observation  about  the  error 
bound  of  the  sparse-dual-based  approach.  The  alternating  iterative 
algorithm  also  converges  fine  for  all  numerical  tests  we  observed. 
The  convergence  guarantee  is  what  we  are  trying  to  establish.  We 
have  made  a  fundamental  observation  that  the  convergence  is  funda¬ 
mentally  a  numerical  stability  issue  of  a  new  “sparse-dual-infused” 
system  matrix  A  (as  in  equation  (2)  in  a  later  description).  The 
next  step  is  to  establish  the  conditions  with  which  “tails”  of  the 
coefficients  are  diminishing. 

[22]  Compressed  sensing  with  frames  is  typically  formulated  to  re¬ 
cover  /  from  the  under-determined  measurements  y  =  Af  (with  or 
without  noise),  assuming  /  =  Dx.  Here  x  is  sparse  and  D  is  a  frame. 
One  typical  and  relatively  successful  (and  natural)  approach  is  to 
simply  write  the  measurement  as  y  =  ADx ,  and  trying  to  recover 
the  coefficient  x  by  various  means.  There  is,  however,  no  good  per¬ 
formance  guarantee  to  such  an  approach,  because  of  the  coupling  of 
A  and  D  in  their  product.  We  have  seen  that  the  sparse-dual-based 
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analysis  approach  can  be  written  as 

(1)  min||n||i  subj.  to  Av  =  b , 

V 

where 

\  AD  1  .  r  I"  b 

(2)  A  =  \l  -  D*(DD*)~lD\  ’  and  h=  [Ax  ' 

Here  Ax  is  the  tail  difference  between  the  canonical  frame  expansion 
coefficients  and  the  sparse  expansion  coefficients. 

We  have  observed  that  such  a  sparse-dual-based  analysis  approach, 
in  the  form  of  (1)  and  (2),  has  a  “decoupling”  functionality  to  decou¬ 
ple  the  product  AD  in  all  traditional  performance  guarantee.  For 
instance,  that  the  unique  i\  solution  is  guaranteed  the  Iq  solution 
as  well  if  and  only  if  A  has  the  Null  Space  Property  (NSP).  NSP  is 
entirely  about  the  null  space  of  A.  In  the  original  problem  of  solving 
x  from  y  =  ADx,  the  NSP  would  be  about  the  null  space  of  AD.  A 
close  examination  shows  that 

ker(A)  =  {v  |  Dv  /  0,  Dv  £  ker(A)}  =  D*(DD*)~1  ( kerA ) . 

As  a  result,  the  kernel  of  A  is  no  longer  about  the  kernel  of  AD,  but 
a  mapping  of  kerA  only.  Hence,  the  sparse-dual-analysis  approach 
decouples  A  from  the  product  with  D  in  the  traditional  approach  of 
solving  for  x  from  y  =  ADx. 

We  believe  that  this  observation  is  capable  to  provide  much  more 
satisfactory  theoretical  insight  about  the  sparse  recovery  problems 
when  signals  are  sparse  with  respect  to  frames. 
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